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^\ • 1 Introduction 

o : 

Problem of cosmological singularity (PCS) is one of the most principal problems of gen- 
eral relativity theory (GR) and relativistic cosmology. The appearance of singular state 
with divergent values of energy density and invariants of curvature tensor in cosmological 
solutions of GR, which limits the existence of the universe in the past, is inevitable accord- 
ing to Hawking-Penrose theorems, if gravitating matter satisfies so-called energy dominance 
conditions [1, 2]. There were many attempts to resolve PCS in the frame of GR as well 
as various generalizations of Einstein's theory of gravitation. The radical idea of quantum 
birth of the Universe was introduced in order to avoid the PCS 3J. A number of regular 
cosmological solutions was obtained in the frame of metric theories of gravitation and also 
other theories, in the frame of which gravitation is described by using more general geometry 
than the Riemannian one (see [4-8] and refs given herein). However, still now we do not 
have the resolution of PCS. Really, the resolution of PCS means not only containing regular 
cosmological solutions, but also excluding singular solutions of cosmological equations by 
using physically reasonable initial conditions. Note that gravitation theory and cosmological 
equations have to satisfy the correspondence principle with Newton's theory of gravitation 
and GR in the case of usual gravitating systems with sufficiently small energy densities and 
weak gravitational fields excluding nonphysical solutions. So metric theories of gravitation 
based on gravitational Lagrangians including terms quadratic in the curvature tensor lead 
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to cosmological equations with high derivatives, and although these theories permit to ob- 
tain regular cosmological solutions with Friedmann and de Sitter asymptotics, however, they 
possess nonphysical solutions also. At last time some regular cosmological solutions were 
found in the frame of superstring theory (brane cosmology), but these solutions also do not 
resolve the PCS (see for example [7-10]). 

Present paper is devoted to study the PCS in the frame of gauge theories of gravitation 
(GTG), at first of all of the Poincare GTG and metric-affine GTG. Note that GTG are nat- 
ural generalization of GR by applying the local gauge invariance principle, which is a base 
of modern theory of fundamental physical interactions (see review [llj). The first attempt 
to apply the simplest Poincare GTG - the Einstein-Cartan theory — in order to resolve 
PCS was made in Refs OEI], where some nonsingular models with spinning matter were 
built. Note that these results depend essentially on classical model of spinning matter, and 
in the case of usual spinless matter the Einstein-Cartan theory is identical to GR and hence 
singular Friedmann models of GR are its exact solutions. The next step to investigate the 
PCS in the frame of Poincare GTG was made in Kel. TTj by using sufficiently general gravi- 
tational Lagrangian including both a scalar curvature and terms quadratic in the curvature 
and torsion tensors. Note that in the frame of GTG quadratic in the curvature tensor terms 
of gravitational Lagrangian do not lead to high derivatives in cosmological equations for ho- 
mogeneous isotropic models modifying Friedmann cosmological equations of GR (see below). 
The conclusion about possible existence of limiting energy (mass) density for usual spinless 
gravitating matter was obtained from deduced cosmological equations leading to regular in 
metrics cosmological solutions. (Later the possible existence of limiting energy density was 
discussed in Ref. jTJ] by modifying cosmological Friedmann equations of GR.) Because cos- 
mological equations of Poincare GTG are valid in the frame of the most general GTG - 
metric-affine GTG I17j. the same conclusions are valid also in metric-affine GTG . Fur- 
ther study of homogeneous isotropic models in GTG has showed that these theories possess 
important regularizing properties and lead to gravitational repulsion effect nearby limiting 
energy density, which can be caused by different physical factors [18-21]. Because the be- 
haviour of cosmological models at the beginning of cosmological expansion in GTG depends 
essentially on equation of state of gravitating matter, we have to know this equation for grav- 
itating matter at extreme conditions in order to build more realistic cosmological models. 
Unified gauge theories of strong and electroweak interactions with spontaneous symmetry 
breaking are the modern base of matter description at extreme conditions. Inflationary cos- 
mology as important part of the theory of early Universe was built by using gauge theories 
of elementary particles [221 EBj- A number of problems of standard Friedmann cosmology 
were resolved in the frame of inflationary cosmology. Most inflationary cosmological models 
discussed in literature are singular and their study is given from Planckian time. Essential 
contribution to energy density in inflationary models is given by so-called gravitating vacuum 
(connected with scalar fields), for which pressure p and energy density p > are connected in 
the following way p = —p. As it was shown in Refs. ^ElEI], in the frame of GTG gravitating 
vacuum with sufficiently large energy density can lead to the vacuum gravitational repulsion 
effect (VGRE) in the case of systems including also usual gravitating matter, that allows 



2 



to build regular inflationary cosmological models. At first the VGRE was discussed in the 
frame of Poincare GTG i n the case of homogeneous isotropic models including radiation 
and gravitating vacuum with p = const. Regularizing role of gravitating vacuum in GTG 
was analyzed in Refs [25-28]. Unlike GR where gravitating vacuum can lead to a bounce 
only in the case of closed cosmological models EDj , in GTG the VGRE allows to build 
bouncing inflationary models of flat, open and closed type. In order to build more realistic 
inflationary models, homogeneous isotropic models including scalar fields and ultrarelativis- 
tic matter were analyzed by applying cosmological equations of GTG |31| . It was shown 
that by certain conditions GTG permit to build regular inflationary cosmological models 
of flat, open and closed type with dominating ultrarelativistic matter at a bounce, and the 
greatest part of cosmological solutions have bouncing character. By taking into account 
that regular in metrics cosmological solutions in GTG take place in the case of gravitating 
matter with p ^ |p, we will study below homogeneous isotropic models including besides 
scalar fields and ultrarelativistic matter also gravitating matter with p ^ |p. In Section 2 
cosmological equations of GTG describing such models are introduced. In Section 3 most 
important general solutions properties of introduced cosmological equations are studied. In 
Section 4 bouncing character of inflationary cosmological models in GTG is analyzed. 



2 Generalized cosmological Friedmann equations in GTG 

Homogeneous isotropic models in GTG are described by the following generalized cosmolog- 
ical Friedmann equations (GCFE) 
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where R(t) is the scale factor of Robertson- Walker metrics, k = +1,0,-1 for closed, flat, 
open models respectively, (3 is indefinite parameter with inverse dimension of energy density, 
M p is Planckian mass, a dot denotes differentiation with respect to time 1 . (The system of 
units with H — c — 1 is used). Eqs.(l)-(2) are identical to Friedmann cosmological equations 
of GR if (3=0. At first the GCFE were deduced in Poincare GTG |14j . and later it was shown 
that Eqs.(l)-(2) take place also in metric-affine GTG [TBIE]. From Eqs. (l)-(2) follows 
the conservation law in usual form 



p + 3H (p + p) = 0, 



(3) 



1 Parameter (3 is defined as (3 = — | (167r) 2 / M 4 , where / is linear combination of coefficients at terms 
of gravitational Lagrangian quadratic in the curvature tensor. 
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where H = j- is the Hubble parameter. Besides cosmological equations (l)-(2) gravitational 
equations of GTG lead to the following relation for torsion function S and nonmetricity 
function Q 

5-ig = ~^ln|l-/3(p-3p)|. (4) 

In Poincare GTG Q = and Eq. (4) determines the torsion function. In metric-affme GTG 
there are three kinds of models in the Riemann-Cartan space-time (Q = 0), in the 
Weyl space-time (S = 0), in the Weyl-Cartan space-time (S ^ 0, Q ^ 0, the function S 
is proportional to the function Q). The value of determines the scale of extremely 

high energy densities. Classical description of gravitational field is valid, if we suppose 
that \(3\~ l < IMp. The GCFE (l)-(2) coincide practically with Friedmann cosmological 
equations of GR if the energy density is small \j3{p— 3p)\ < 1 (p / The difference 

between GR and GTG can be essential at extremely high energy densities \/3(p — 3p)| > 1. 
Ultrarelativistic matter with equation of state p — |p is exceptional system because Eqs. 
(l)-(2) are identical to Friedmann cosmological equations of GR in this case independently 
on values of energy density. The behaviour of solutions of Eqs. (l)-(2) depends essentially on 
equation of state of gravitating matter at extreme conditions and on sign of parameter /3. By 
given equation of state p = p(p) the GCFE (l)-(2) lead to regular in metrics cosmological 
solutions in the following cases: 1) > and p < hp, 2) (3 < and p > |p [I]. The 
investigation of models including scalar fields on the base of Eqs (l)-(2) shows that the 
choice (3 < permits to exclude the divergence of time derivative of scalar fields |ITj. In 
connection with this we put below that parameter f3 is neg ative and \f3\- 1 < 1 • Af*. 

In the case of gravitating vacuum with constant energy density p v = const > the 
GCFE (l)-(2) are reduced to Friedmann cosmological equations of GR and S — Q — 0, this 
means that de Sitter solutions for metrics with vanishing torsion and nonmetricity are exact 
solutions of GTG and hence inflationary models can be built in the frame of GTG. 

In order to analyze inflationary cosmological models in GTG let us consider systems 
including scalar field minimally coupled with gravitation and gravitating matter in the form 
of several components with energy densities pi and pressures pi = WiPi, where Wi = const (i 
is a number of component). Later we will consider two components: ultrarelativistic matter 
with energy density p r and w r = |, and gravitating matter with w > |. If the interaction 
between scalar field and all components of gravitating matter is negligible, the energy density 
p and pressure p take the form 

P = \fi + V{<j>) + P = - V(<f>) + J2 w iP^ ( 5 ) 

where the summation is taken over all components. The conservation law (3) leads to the 
scalar field equation 

dV 

<P + 3H<P=- — , (6) 

where V is scalar field potential and to the conservation laws for all components of gravitating 
matter 

p\ + 3H Pi (l + Wi ) = 0. (7) 
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Eqs (7) have integrals in the following form 

PtR 3(l+m) = congt 

By using Eqs.(5)-(7) the GCFE (l)-(2) can be transformed to the following form 
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Relation (4) takes the 



Unlike GR the cosmological equation (9) leads to essential restrictions on admissible values 
of scalar field and gravitating matter. Imposing (3 < 0, we obtain from Eq. (9) in the case 
k = 0,+l 

Z>0 or ft <AV +\(3\~ 1 -^piiZwi-l) . (12) 

Inequality (12) is valid also for open models discussed below. The region E of admissible 
values of scalar field 0, time derivative <fi an d energy densities pi (excluding p r ) in space P 
of these variables determined by (12) is limited by bounds L± 



= ± 



From Eq. (9) the Hubble parameter on the bounds L± is equal to 



(13) 



According to (14) the right-hand part of Eq. (11) is equal to |if , this means that the torsion 
(nonmetricity) will be regular, if the Hubble parameter is regular. In Sections 3-4 our main 
attention will be turned to study properties of solutions of GCFE (9)-(10). 



3 Some general properties of GCFE solutions for in- 
flationary models 

Let us consider the most important general properties of cosmological solutions of GCFE 
(9)-(10). At first, note by given initial conditions for scalar field ((f), <p) and values of R 
and pi there are two different solutions corresponding to two values of the Hubble parameter 
following from Eq. (9): 

H = W> ± VP 

± l-2/?(2V + ^)-§0$>(9u;?-l)' 

where 

^=^{E^ + ^ 2 + ^-^K-^-Ea(3^-1)] 2 }^-^^ 2 >0 (16) 

Unlike GR, the values of H + and H- in GTG are sign-variable and, hence, both solutions 
corresponding to H + and H_ can describe the expansion as well as the compression in 
dependence on their sign. Below we will call solutions of GCFE corresponding to H + and 
H_ as i/ + -solutions and iJ„-solutions respectively. In points of bounds L± we have D — 0, 
H + = H- and the Hubble parameter is determined by (14). Eqs. (9)-(10) are satisfied on the 
bounds L±, corresponding solutions of GCFE - L-t-solutions - are their particular solutions; 
i/_-solutions reach the bounds L± and //^.-solutions originate from them (see below). By 
using Eqs. (15)— (16) it is easy to shown, that in points of bounds L± the derivatives H + 
and H_ are equal: \imH + = lim.ff_ at Z — > 0. As result, we have the smooth transition 
from i/_-solution to i7 + -solution on bounds L±. At the same time, the value of the time 
derivative of H for L^-solutions according to (14) is not equal to limH± at Z — > 0, and by 
transition from if_-solution to L±-solutions and from L±-solution to iJ + -solution a finite 
jump of the derivative H takes place. Note, that according to Eq. (11) the functions S and 
Q have the following asymptotics at Z — > 0: \S — \Q\ ~ Z^ 1 / 2 for H + - and if_-solutions. 

In order to study the behaviour of cosmological models at the beginning of cosmo- 
logical expansion, let us analyze extreme points for the scale factor R(t): Rq = R(0), 
H = H(0) = 0. Denoting values of quantities at t = by means of index "0", we obtain 
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from (9)-(10): 
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the value of Hq is positive. By using Eq.(17) we can rewrite the expression of Hq in the form 



#o 



7T 



Vn 



'1 - w%) 



4 



p i0 {3wi - 1) 



+3f3(V" <Pl - V 



r2\ 



2k 



}{l-2/? 2Vb + ^ + ^p iQ (9^ 2 -l) } 



(19) 



We see from (19) unlike GR the presence of gravitating matter (with Wi < 1) does not 
prevent from the bounce realization 2 . In the case of various scalar field potentials applying in 
inflationary cosmology Eq.(17) determines in space P so-called "bounce surfaces" depending 
on parameter j3 and energy density of ultrarelativistic matter parametrically. In the case of 
closed and open models families of bounce surfaces depend also on the scale factor i?o-By 
giving concrete form of potential V and choosing values of R , <p , (j) and p i0 at a bounce, we 
can obtain numerically particular bouncing solutions of GCFE for various values of parameter 
(3. 

The analysis of GCFE shows, that properties of cosmological solutions depend essentially 
on parameter (3, i.e. on the scale of extremely high energy densities. From physical point of 
view interesting results can be obtained, if the value of \(3\~ l is much less than the Planckian 
energy density [33], i-e. in the case of large in module values of parameter (3 (by imposing 
M p = 1). In order to investigate cosmological solutions at the beginning of cosmological 
expansion in this case, let us consider the GCFE by supposing that 
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2 In GR a bounce is possible only in closed models if the following condition Vq — 4% — \ EJ Pio(l + 3u>i) > 
takes place. 
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Note that the second condition (20) does not exclude that ultrarelativistic matter energy 
density can dominate at a bounce. We obtain: 
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Eqs. (21)-(22) do not include radiation energy density, which does not have influence on the 
dynamics of inflationary models in the case under consideration (although, as it was noted 
above, the contribution of ultrarelativistic matter to energy density can be essentially greater 
in comparison with scalar field and other components of gravitating matter), moreover Eqs. 
(21)-(22) do not contain the parameter f3. According to Eq. (21) the Hubble parameter in 
considered approximation is equal to 
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and extreme points of the scale factor are determined by the following condition 
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From Eq. (22) the time derivative of the Hubble parameter at extreme points is 
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or according to Eq. (24) we can rewrite the expression Hq in the following form 
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Obviously Eqs. (24)-(25) correspond to (17)-(18) in considered approximation. 

Now in order to investigate inflationary cosmological models in GTG we will analyze 
models including scalar fields with positive potentials 3 V > 0, ultrarelativistic matter and 
component of gravitating matter with energy density p\ and w\> \. 



4 Analysis of inflationary cosmological models in GTG 

At first let us consider models including scalar field and ultrarelativistic matter studied in 
Ref. [HI]- In this case the space P is reduced to the plane of variables (0, 0), bounds L± — to 

two curves <fi = ± (4V + |/9| -1 ) 2 and bounce surfaces — to corresponding bounce curves on 
this plane. As it was shown in Ref. [3T] the greatest part of inflationary cosmological solutions 
in the case under consideration have regular bouncing character, although singular solutions 
exist also because of divergence of particular L±-solutions. To analyze regular bouncing 
solutions we have to examine bounce curves and ff±-functions defined by Eq.(15) on the 
plane P. Bounce curves have simple form, if the scale of extremely high energy densities 
is much smaller than the Planckian energy density. Then according to (24) bounce curves 
do not depend on parameter (3 and bounce curves of flat models are two curves B\ and B 2 
determined by equation 

on the plane P, which are situated near bounds L + and L_ respectively. 4 Each of two 
curves B 12 contains two parts corresponding to vanishing of H + or H_ and denoting by 

3 We do not consider in present paper negative scalar field potentials [55] , 

4 The neighbourhood of origin of coordinates on the plane P is not considered in this approximation, the 
behavior of bounce curves near origin of coordinates was examined in Ref. |84| , where scalar fields superdense 
gravitating systems were discussed. 
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(Bi + , B 2+ ) and (B^, B 2 -) respectively. If V is positive (negative) in quadrants 1 and 4 
(2 and 3) on the plane P, the bounce will take place in points of bounce curves B x+ and 
B 2 + {Bi- and B 2 J) in quadrants 1 and 3 (2 and 4) for i/ + -solutions (//--solutions) (see 
Fig. 1). To analyze flat bouncing models we have to take into account that besides regions 




Figure 1: Bounce curves for i/+-solutions and //--solutions in the case of potential V = 

lying between bounds L± and corresponding bounce curves on the plane P the sign of values 
H + and H for applying potentials is normal: H + > 0, H < 0. The Hubble parameter H + 
is negative in regions between curves (L+ and Bi + ), (L_ and B 2+ ), and the value of i/_ is 
positive in regions between curves (L + and -Bi_), (L_ and B 2 ). Any regular cosmological 
solution has to include both i/_- and // + -solution. The regular transition from //^-solution 
to i/ + -solution takes place on the bounds L± where H + = //_. If H + - and //_-solution have 
with L_|--bounds one common point, corresponding bouncing solution is regular in metrics, 
the Hubble parameter H and its time derivative. If H + - and //--solution are glued with 
particular L + - or L_-solution in different points of bounce, corresponding bouncing solution 
is regular in R{t) and H{t), but the derivative H in points of gluing has a finite jump. In 
the case of open and closed models Eq. (24) (pi = 0, with the exception p r 7^ 0) determines 
1-parametric family of bounce curves with parameter R . Bounce curves of closed models 
are situated on the plane P in region between two bounce curves Bi and B 2 of flat models, 
and in the case of open models bounce curves are situated in two regions between the curves: 
L + and B\, L_ and B 2 . Because the behaviour of bounce curves for open models is like to 
that for flat models, the situation concerning bouncing inflationary solutions in the case of 
open models is the same as described above situation for flat models. Unlike flat and open 
models, for which H + = H_ only in points of bounds L± and regular inflationary models can 
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be built if H + — and //.-solutions reach bounds L±, in the case of closed models the regular 
transition from //.-solution to // + -solution is possible without reaching the bounds L±. It 
is because by certain value of R according to (16) we have H + = //_ in the case, if Z ^ 0. 
Such models are regular also in torsion and/or nonmetricity. Regular inflationary solution 
of such type was considered in Ref.jHBj- 

However, in discussed case there are also singular solutions. At first note that particular 
L±-solutions are singular. Because bounds L± for applying scalar field potentials tend to 
infinity, we have that scalar field satisfying on the bounds the equation <fi = 2V diverges 
in the past (quadrants 2 and 4) and in the future (quadrants 1 and 3). In consequence 
of this any solution including /^-solution (or //.-solution) glued with one of the bounds 
L± is singular. Hence the problem of excluding of singular solutions is connected with 
regularization of particular L±-solutions. 

In general case, when approximation (20) is not valid, bounce curves of cosmological 
models including scalar fields and ultrarelativistic matter determined by Eq.(17) depend 
on parameter (3 . By certain value of f3 we have 1-parametric family of bounce curves 
with parameter p r0 for flat models, and we have 2-parametric families of bounce curves for 
closed and open models with parameters Rq and p r Q. The situation concerning cosmological 
solutions of Eqs. (9)-(10) does not change. 

Now let us analyze models including gravitating matter with energy density p\ and 
W\ > | besides scalar field and ultrarelativistic matter. The bounds L± in this case depend 
on energy density p\ and they change by evolution of models 



= ±yj4V + \f3\~ 1 - Pl (3wi - 1) . (27) 
The Hubble parameter on bounds L± according to (14) is equal to 

H = W (2g) 

I - 2(3(2V + ft) - lp Pl {9wl- 1)' 

Then the equation of scalar field (6) on bounds L± because of (27)— (28) takes the following 
form 

1 - 48V - 

If | < w± < 1 and (3 < 0, the denominator of a fraction in (29) is positive, and for various 
potentials V applying in inflationary cosmology (in particular, V = ^m 2 (j) 2 , V = \\(j) 4 
etc.) Eq. (29) describes finite variations of scalar field between positive maximum and 
negative minimum values of <fi at <p = (see Fig. 2). This means that particular solution of 
GCFE becomes regular, and it has bouncing character; a bounce takes place in points where 
(p = (according to (14) and (19) we have in these points H = and H > 0). As result, 
all inflationary cosmological solutions of GCFE are regular. The analysis of inflationary 
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solutions reaching the bounds L± by numerical integration of Eqs. (11) and (6) is difficult, 
because the coefficient at H in Eq. (11) tends to zero at Z — > 0. Note that bouncing 
character have solutions not only in classical region, where scalar field potential, kinetic 
energy density of scalar field and energy density of gravitating matter do not exceed the 
Planckian energy density, but also in regions, where classical restrictions are not fulfilled 
and according to accepted opinion quantum gravitational effects can be essential. 
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Figure 2: Numerical solution of Eq. (29) for V = \m 2 (j) 2 
(m = 10" 6 M p , w = l,P= -10 6 M- 4 ). 

As illustration of obtained results we will consider particular bouncing cosmological in- 
flationary solution for flat model by using scalar field potential in the form V = \m 2 4> 2 
(m = 10~ 6 M P ). The solution was obtained by numerical integration of Eqs. (6), (10) and by 
choosing in accordance with Eq.(17) (or (24)) the following initial conditions at a bounce: 
0o = v^lO 3 M p , 0o = 10~ 3 M 2 , p w = 1.4999 • 10" 6 M 4 (w = 1, (5 = -10 14 M" 4 ); in accor- 

R 6 

dance with (8) the formula Pi(t) = Pio^^y was used, initial value of i? can be arbitrary. 

A bouncing solution includes: quasi-de-Sitter stage of compression, the stage of transition 
from compression to expansion, quasi-de-Sitter inflationary stage, stage after inflation. The 
dynamics of the Hubble parameter and scalar field is presented for different stages of ob- 
tained bouncing solution in Figures 3-5 (by choosing M p = 1). The transition stage from 
compression to expansion (Fig. 3) is essentially asymmetric with respect to the point t — 
because of 0o 7^ 0. In course of transition stage the Hubble parameter changes from maxi- 
mum in module negative value at the end of compression stage to maximum positive value 
at the beginning of expansion stage. The scalar field changes linearly at transition stage, 
the derivative grows at first from positive value 0i ~ 1.6 • 10~ 7 to maximum value ~ O 
and then decreases to negative value 2 ~ — 1.6 • 10~ 7 . Quasi-de-Sitter inflationary stage and 
quasi-de-Sitter compression stage are presented in Fig. 4. Although the GCFE (9)-(10) and 
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Figure 3: The stage of transition from compression to expansion. 
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Figure 4: Quasi-de-Sitter stage of compression and inflationary stage. 



their approximation (21)-(22) have different structure from cosmological Friedmann equa- 
tions of GR, like GR the time dependence of functions H(t) and <p(t) at compression and 
inflationary stages is linear. The amplitude and frequency of oscillating scalar field after in- 
flation (Fig. 5) are different than that of GR, this means that approximation of small energy 



densities 



(AV - 2 - 2 Pl ) 



C 1 at the beginning of this stage is not valid; however, the 



approximation (21)-(22) is not valid also because of dependence on parameter (3 of oscilla- 



13 



tions characteristics, namely, amplitude and frequency of scalar field oscillations decrease by 
increasing of |/3| [33]. The behaviour of the Hubble parameter after inflation is also nonein- 
steinian, at first the Hubble parameter oscillates near the value H = 0, and later the Hubble 
parameter becomes positive and decreases with the time like in GR. Before quasi-de- Sitter 
compression stage there are also oscillations of the Hubble parameter and scalar field not 
presented in Figures 3-5. Ultrarelativistic matter, which could dominate at a bounce as 
well as another component of gravitating matter have negligibly small energy densities at 
quasi-de Sitter stages. At the same time the gravitating matter could be at compression 
stage in more realistic bouncing models, and scalar fields could appear only at certain stage 
of cosmological compression. As it follows from our consideration regular character of such 
inflationary cosmological models has to be ensured by cosmological equations of GTG. 



10 7 H 10 2 (J) 




Figure 5: The stage after inflation. 



5 Conclusion 

As it is shown in our paper, GTG permit to build regular inflationary cosmology, if gravitat- 
ing matter at extreme conditions satisfies the following restriction |p < p < p and indefinite 
parameter j3 is negative. All inflationary cosmological solutions for flat, open and closed 
models are regular in metrics and the Hubble parameter. The presence of scalar fields lead- 
ing to appearance of inflation in cosmological models changes essentially the structure of 
GCFE, as result a family of closed models regular in metrics, the Hubble parameter and 
torsion and/or nonmetricity appears. To build realistic inflationary cosmological models not 
limited in the time we have to know the change of equation of state of gravitating matter 
by evolution of the Universe. 
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